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Abstract 

Quantum 3D R - matrix in the classical (i. e. functional) limit gives a symplectic map of 
dynamical variables. The corresponding 3-D evolution model is considered. An auxiliary problem 
for it is a system of linear equations playing the role of the monodromy matrix in 2D models. 
A generating function for the integrals of motion is constructed as a determinant of the auxiliary 
system. 
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1 Introduction 



In this paper we formulate a classical discrete space - time model (DSTM) 1 , which is in some sense the 
classical counterpart of the Zamolodchikov - Bazhanov - Baxter 3D lattice spin model |I| || |j. One 
way to obtain a classical model associated with ZBB one is to consider the path integration motivated 
quasi-classical limit of ZBB, this was done in Ref. Q. 

But here we choose another way: we consider a symplectic map given by the functional limit of the 
operator 3D R - matrix ||, In this approach we elude the usage of the Hirota - Miwa equation, 
because of this map is the local symplectic map for three pairs of dynamical variables. This map we use 
to construct 31? DSTM with dynamical variables belonging to the vertices of the kagome lattice. 

An auxiliary linear problem for such maps in the operator as well as in the functional approaches 
were found recently in fl~o|| . This problem is an over - defined set of linear equations with compatibility 
conditions giving the map. Analogous system of linear equations for whole two dimensional lattice 
plays the role of the monodromy matrix in 2D models, and remarkably, the determinant of this system 
corresponds to the trace of the monodromy matrix and gives a generating function for integrals of motion. 

In the main part of this paper (section 2: "The model") we formulate the evolution model, introduce 
linear variables and show how the integrals of motion can be derived. An interpretation of the results 
and some aspects of further investigations are given in section 3: "Discussion". 

Section 2 contains all the essential results of the paper, but in such short form they look artificially. 
All cumbersome details concerning the derivation of the model, its geometrical interpretation and correct 
application of the original linear problem from Ref. |l0[ are given in Appendix. 

2 The model 

2.1 Algebraic definition of the evolution 

Consider a Z 2 span of a two dimensional plane generated by two independent shifts a and b such that 
point (n,m) of the span we present as P — a n ■ b rn . Assign to each point P three pairs of dynamical 
variables: 

[u x {P) : wx{P), u 2 (P),w 2 (P), u 3 (P),w 3 (P)\. (2.1) 
For this discrete space system introduce a discrete time evolution 

[u k (n,P),w k {n,P)] ' ^ [u k {n + l,P),w k {n + l,P)] , k -1,2, 3, (2.2) 

defined by the following relations: 

wi(n, P) w 2 {n, P) + w 2 (n, P) u 3 (n, P) + k 2 u 3 (n, P) w 3 (n, P) 

w 3 (n,P) 

(2.3) 

P) u 2 (n,P) w 2 (n,P) 
ui(n, P) w 2 (n, P) + u 2 (n, P) w 3 (n, P) + k 2 ui(n, P) w 3 (n, P) 



Rcfs. [hi ti. hi arc the rather uncompleted list of the publications which are to be mentioned 



Wi(n+l,P) = 

< 

ui(n+l,P) = 



1 



W2(n + 1, a ■ P) 
112(11 + 1, a ■ P) 

w 3 (n + l,b-P) 
u 3 (n + 1,6-P) 



w\(n,P) 1112(11, P) w 3 (n,P) 



wi(n, P) 11)2(11, P) + w 2 (n,P) 113(11, P) + n 2 113(71, P) w 3 (n,P) 

u±(n,P) u 2 (n,P) u 3 (n,P) 

Wi(n, P) u 2 (n, P) + u 2 (n, P) u 3 (n, P) + n 2 Ux(n, P) Wi(n, P) 

u 2 (n,P) w 2 (n,P) w 3 (n,P) 

Ui(n, P) w 2 (n, P) + u 2 (n, P) W3(n, P) + k 2 U\(n, P) w 3 (n, P) ' 

w\(n,P) u 2 (n,P) + u 2 (n,P) u 3 (n,P) + n 2 u\(n, P) w\(n, P) 

m(n,P) 



(2.4) 



(2.5) 



Here k 2 is an arbitrary parameter common for all P. The main property of this evolution is that it 
conserves the Poisson brackets 

{u k (n, P) , w m (n,Q)} = u k (n,P) w m (n,Q) 6k, m &P,Q , (2-6) 



i. e. map (2.2) is the symplectic map. 



2.2 Auxiliary linear problem 

Introduce a set of auxiliary variables ipl(n, P) assigned to the points of the 2D lattice. Let the evolution 

otiP* k {n,P)be 

iPt(n + l,P) = r 3 (n,P), r 3 (n + l,P) = ^(n, a^b' 1 ■ P) , (2.7) 



and 



r 2 (n + l,P) = ri (n,P)n( - ^ ) 

Ui(n,P)w 2 (n,P) u 2 (n, P) w 3 {n, P) 

-r 3 (n,a-P) 

+r 3 (n,p) 



Ui(n + l,P) w 3 (n,P) 
113(11, P) 



U2(n + 1, a ■ P) 



Thus i/jt(n,P) are the linear variables. Define functions J£(n,P) as follows: 

J*(n, P) = ipl(n, P)wt(n, P) ~ ip^(n, P)ui(n, P) + ip^(n, aP)n + $J(n, bP)nui(n, P)wi(n, P), 

J 2 *(n, P) = iPl(n, P)k + ip* 2 (n, aT x P)KU2{n, P)w 2 (n, P) + $}(n, P)w 2 (n, P) - ip^(n, bP)u 2 (n, P), 

J 3 *(n, P) = -i\>l(n, P)u 3 (n, P) + i/j* 2 (n, b- 1 P)w 3 (n, P) + iP* 3 (n, aP) K + iP* 3 {n, P)nu 3 (n, P)w 3 (n, P). 

(2.9) 
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One can verify easily 



^ iti (n + 1, PWi (n + 1, P) , k 

J x (n + 1, P) = J 2 * (n, P) TzhL. -7Z o^ + J 3 (*> p ) 



u 2 (n,P)w 2 (n,P) ' w 3 (n,P) 

j;(n + l,a.P) = j; { n,P) ^1 + J?(n,P) , 

ui(n, P) w 3 (n,P) 

r 3 (n + l, b .P) = -JKn,P) + J 2 >,P) «a(" + 1> PW» + 1,6 • P) 

ui(n,P) U2{n, P)W2(n, P) 

(2.10) 

Regarding ipt(n, P) and Jt(n, P) as components of row vectors and J*(n), we obtain three formal 

matrix equations: first, Eq. (|2.9|) we rewrite as 



J*(n + 1) = $*(n) • D(n) , (2.11) 

second, Eqs. fl2.7||2.g| ) we rewrite as 

**(n + l) = **(n) ■ K(n) , (2.12) 



and finally, Eq. (2.1C) - as 

J*(»+l) = J*(n) • M(ra) . (2.13) 

Hence 

J*(n + 1) = #*(n) • K(n) • D(n + 1) = **(n) • D(n) • M(n) , (2.14) 
so the main relation of this subsection arises: 

K(n) ■ D(n+1) = D(n) • M(n) . (2.15) 



Note, there is some ambiguity in the definition of K and M. Namely, for (2.7) given ip%(n + 1,P) 
can be obtained in general from 

Jf(n+l,P) = Ji(n,P) a(n,P) + J 2 * (n, P) /3(ra, P) + J 3 *(n, P) 7 (n, P) (2.16) 

with arbitrary a, /?, 7. For K given - M can be restored directly. 

2.3 Finite lattice and integrals of motion 

Reduce now the infinite Z 2 span of the two dimensional plane to the finite 7? N span of a two dimensional 
torus. Namely, let 

u k (n,a N -P) = u k (n,b N -P) = u k {n,P), 

w k (n,a N -P) = w k (n,b N -P) = iu fc (n,P), (2.17) 

k = 1,2,3, for all P. Obviously, these conditions are conserved by evolution (^3-2.5). Pure periodic 
boundary conditions are natural for the dynamical variables, while for auxiliary ip k (n, P) one can impose 
quasiperiodical boundary conditions: 

Win, a N -P) = r k (n, P) x , $S(», b N ■ P) = ffi(n, P) y . (2.18) 
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Variables x, y play the role of the spectral parameters. Now all matrices K,M,D arc finite (3A 2 x 37V 2 ) 
dimensional matrix functions on x,y. 



Take now the determinants of the left and right hand sides of Eq. (2.15). It is easy to see that 

dct K(n, X, y) = x- N y- N det K(n, 1, 1) = x~ N y~ N K{n) , (2.19) 

and 

detM(n,x,y) = x~ N y~ N det M(n, 1, 1) = x~ N y~ N M(n) . (2.20) 

Denote further 

D(n,x,y) = detT>(n,x,y) . (2.21) 



Hence Eq. ( 2.\l ) gives 



D(n + l,x,y) = ^^-D{n,x,y). (2.22) 
Taking the series decomposition of D, 

D(n,x,y) = ]T x l y m D l<m (n) , (2.23) 

l.nt 

we obtain the integrals of motion as the ratios of -Dfc,m* 

Ii, m [u k (n,P),w k (n,P)} = n A,m( "\ , (2.24) 

where li^ m are functionals of Mfe(n, P), Wfe(n, P) and 

I l<m [u k (n,P),w k (n,P)] = I hm [u k (n+l,P),w k {n + l,P)}. (2.25) 

3 Discussion 

So, we have described the completely integrable classical discrete space - time model. This model can 
be quantized (see Appendix for the details) , and general approach giving relation ( |2.15 ) has a quantum 



counterpart. This gives a way to construct a 3D analogue of the Bethe ansatz. Several papers concerning 
eigenvectors and eigenfunctions for 3D transfer matrices and evolution operators are to be mentioned 
here : and p|. 

Mention also the evolution operator for the affine Toda field theory (ATFT), given in Ref. [fl6|| . 
First, the classical limit of ATFT is closely connected with our evolution: the limit case when k 2 = 

and —^t — }■ — — t — — — - i — ► for all P. turning the frame of reference n, a, b appropriately and hiding 
Ui{P) wi(P) w 2 {P) 

one space direction into the rank of one can obtain a map corresponding to ATFT. Second, the 

approach to 3D models as to 2D models with high rank symmetry group implies the tedious technique 
of the nested Bethe ansatz while we try to avoid this using "relativistic" approach. 

2 This operation in terms of classical variables is defined badly, and in operator language gives some central elements. 
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Appendix 

A Definition of map R 

Consider the following map: 

[ui,wi;u 2 ,w 2 ;u 3 , w 3 ] t 
where Uk, Wk are the following functions of Uk,Wk'- 



[Ui,Wi 1 U2,W 2 ;U 3 ,W 3 \ 



(A.l) 



Wi 
11)2 

w 3 

W2 

w 3 



Wi w 2 - 


\- W 2 U 3 + 


K 2 


u 3 w 3 


w 3 

W± W 2 W 3 


W\ w 2 - 


f- W2 U 3 + 


K~ 


u 3 w 3 




1l 2 W 2 W 3 






U\ W 2 " 


f- U2 W 3 + 


K 2 


U\ w 3 


is given by 








ih\ w 2 u 3 






Wi w 3 - 


h u 3 w 3 + 


K- 


Wi w 2 


W\ w 3 - 


h U 3 W 3 + 


K : 2 


Wi w 2 


ill W2 " 


u 3 

f U 2 U>3 + 




U 2 W2 




Ul 







Ul 

u 3 

Ul 
»2 
"3 



Ui U 2 W 2 



Ul U)2 - 


f u 2 w 3 + 


K 2 Ul W 3 




Ul U 2 U 3 




Wi u 2 - 


f M2M3 + 


K 2 Ul Wi 


Wi U2 ~ 


f U2U 3 + 


K 2 Ul Wi 




Ul 




Ul W2 


f u 2 w 3 + 


K 2 U 2 U>2 




w 3 




Ul Wi 


-\- ill u 3 + 


K 2 ll 2 ll 3 




Wi 






Wi u 2 u 3 




Ul Wi 


f M1U3 + 


K 2 ll 2 ll 3 



(A.2) 



(A.3) 



The reader can easily recognize R in our evolution formulae (2^-2^5) up to some re-enumeration of final 
states' coordinates. The map R can be extracted from the operator function R (see ||). In few words, 
in the operator formulation R is a function of three independent commutative pairs u^, w^, k = 1, 2, 3, 
such that 

Ufe • w fc = q w fe • u fe , (A. 4) 

and 

ii fe = R u fc R 1 , w fc = R w fe • R- 1 . (A.5) 

Map ( |A.2j ) appears in the limit q — y 1, q 1 ^ 2 — > — 1, u — > u, w — > w. As the consequence of ( |A.4[ ) and 
( |A.5| ), R conserves the Poisson brackets (pT^). 

Note that the map R obey the Functional Tetrahedron Equation |l^, |l^, |1| , this fact provides 
the integrability. 

Few words concerning the structure of R. Introduce useful variables u,w, u, w, s as follows: 



_ u Ui 

w = Wi W2 , u — U2 u 3 , s — = — 

w w 3 



Then R conserves u, w and s: 



w 3 Wi 
u = — , w = — 
w 2 u 3 



u = u , w = w 



U Ui 

w u 2 



(A.6) 
(A.7) 

(A.8) 



G 



and R acts on u, w as follows: 

I 1 + w + k 2 u , s 1 + w + k 2 u , . « 

u = — , w — — 5— . (A. 9) 

U S + W + S K z U W 1 + W + S K z U 

Define the function H: 

H(u,w) = w + — + k 2 (1 + w) (- + u — ) . (A. 10) 

w u w 

It can be easily verified, 

H(u,w) = H(u,w) . (A.ll) 

Note, H (u, w) = h is a genus one elliptic curve. Thus H (u, w) can be considered as the Hamiltonian, 
corresponding to R, and the evolution 

H(u,w) r^ili = {H(u,w), u} , H(u, w) T ^ll = { H (u, w) , w } (A.12) 

AT CLT 



gives map ( A.9) as the first rational point. 



B Evolution 



Define now the geometrical evolution model with the help of the map (A.1,A.2). 

Consider the kagome lattice formed by 3A^ 2 oriented lines as it is shown in Fig. A. 




Fig. A. A fragment of the kagome lattice. 



The lattice consists of two types of triangles and one type of hexagons. Choose the triangles marked in 
Fig. A for the enumeration of the lattice elements. A shift of the triangle P to any other triangle of 
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the same type is s superposition of two independent shifts, and because of the existence of three similar 
directions of the kagome lattice it is convenient to use the homogeneous notation for a multiplicative 
shift operator 

7f -4 ■ tZ : a + /3 + 7 = 0, (B.l) 

with the action of t^ x t 2: t{~ T3 and t% t$ defined in Fig. A. 

Consider now triangle P in details. Enumerate the vertices of P by the numbers 1, 2 and 3. Define 
the evolution U of the lattice as the mutual change of all the P - type triangles in the spirit of the Yang 
- Baxter equivalence as it is shown in Fig. B. 




Fig. B. The action of U. 



The lattice obtained is again the kagome lattice with the similar set of P - type triangles marked as tu P 
in Fig. B. We regard this change as the evolution from time n to time n + 1, and obviously the time 
can be extracted from P as a power of Tk (recall, translations (B.l) for a time fixed are of zero power). 
Thus three shifts t±,t 2 ,t 3 are the elementary shifts of a three dimensional rectangular lattice, and 
the evolution U is the saw transfer matrix of the 3D lattice (see Refs. 14, fL5| for geometrical details). 



Assign now the pair of dynamical variables Uk(P),Wk(P) to vertex k of triangle P. 
evolution for these variables via the map R jATpH ) as follows: 



U 



R 



[u 1 (P),w 1 (P) ] u 2 (P),w 2 {Py, u 3 (P),w 3 (P)} 
[ui (n P) , wi (ti P) ; u 2 (r 2 P) , w 2 (r 2 P) ; u 3 (r 3 F) , w 3 (r 3 P)] 

Denote now for the shortness the translations for the time fixed as 

-i _ b = _-i _ 



Define the 



(B.2) 



(B.3) 



« — '1 '2 ) " — '1 '3 ) 

and exhibit the "time" as a power of t\ in the decomposition of three dimensional P with respect to 
Ti,a,b: 

P = rf • a a ■ b 13 ^ (n, P = a a -b 13 ) . (B.4) 



with this definition of time and space coordinates, map (B.2) gives exactly the evolution (2.3-2.5). 
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C Basic auxiliary linear problem 



A version of an auxiliary linear problem for 3D models was suggested in Rcf. [Q. Here we recall it 
briefly for the case of commutative dynamical variables and show how integrals of the motion can be 
derived with a help of it. 
We start from the basic definition. 

Consider some vertex V from our (kagome) lattice. Complete pair u, w assigned to V by an "internal 
current" while observable (or "external") currents we are intending to measure in the faces. Internal 
current as produces additive contributions into four external currents of four faces surrounding the 
vertex. The values of the contributions are shown in Fig. C. 




Fig. C. The current vertex. 



The rules of the game are quite simple: 

• any observable face current is the sum of the surrounding vertices contributions, and 

• the face current for any closed face is zero. 

This indeed is the decent linear problem due to the following remarkable feature: for a given geometry 
of a graph the relations between the outer face currents (with the conservation laws being taken into 
account) define the dynamical variables of the elements of the graph unambiguously. 

Consider now the problem of the equivalence of two Yang - Baxter - type graph as in Fig. B. In Fig. 
D these graphs are redrawn with the face currents' contributions marked out. 
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Fig. D. The equivalence. 



First, the conservation laws for the closed faces <f> a = <j>h = 0. Second, the equivalence means that the 
outer currents <pb, <f>g for the left ang right hand side graphs coincide. Thus restoring the values of the 
vertex contributions into the face currents, one obtains the following system: 

(f) h := wi 4>\ + k cf> 2 - u 3 (f> 3 = , 

4>a ■= - U\ <j)i + k u 2 w 2 4>2 + w 3 4> 3 — , 



(f>b ■■= wi 4>i — w 2 4>2 + k u 3 w 3 (f> 3 , 

4> c ■= k fa = k 4>i + k 4> 3 , 

4>d ■= - u 3 03 = k ui m (f>i - u 2 (f>2 , 

4>e ■= - Ul (f>i = - U2 (j>2 + k (f> 3 , 

(f>f : = k u 2 w 2 <p2 — k til w\ <pi + k u 3 w 3 (f> 3 , 
4> g := w 3 4> 3 — k cpi + w 2 4>2 ■ 

(C.l) 

This system is eight homogeneous linear equations for six internal currents 4>k,4>k- The initial graph 
(let hand side of Fig. D) is described by two independent currents, so the rank of 8 x 6 matrix of the 
coefficients must be 6 — 3 = 4. This is the zero curvature condition, and the solution of it is given by 

w\ = w 2 A 3 , tii = A^ 1 w^ 1 , 
w 2 = Ag 1 w 1 , u 2 = A^ 1 u 3 , 

w 3 = A^ 1 u^ 1 , u 3 = u 2 A x , (C.2) 
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with 

Ai = u^ 1 u 3 + u^ 1 w 1 + k w 1 v>2 ^ , 

A 2 = u^ 1 w^ 1 + u^ 1 ■ w^ 1 + k 2 u^ 1 w^ 1 , 

A 3 = w x Wg 1 + w^ 1 u 3 + n 2 w^ 1 u 3 . (C3) 

This is exactly the map R. 

D Linear problem for the whole lattice 

Consider now the current system for the whole kagome lattice. All faces are closed, so all the face 
currents are zeros. 

Describe the system of "vacuum currents" in details. As it was mentioned, the lattice consists of two 
types of the triangles and one type hexagons. The elementary "puzzle" is shown in Fig. E. 



V 2 {aP) V 3 (bP) 




Fig. E. An elementary puzzle. 
Thus for each "puzzle" three "vacuum" equations are: 

MP) = f Wl (p)MP) + *MP) - MP) MP) = o (D.i) 

MP) d = -MP) MP) + KU 2 (aP)w 2 (aP) MaP) + w 3 {bP) fo(bP) = (D.2) 

MP) d = w 2 (p)Mp) + kMp)Mp)Mp) + kMMp) + k mMp) 
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+KUi(b- 1 P)wi(b- 1 P)(f>i(b- 1 P) -u 2 {b- 1 P)(f)2(b~ 1 P) = 



(D.3) 



We are looking for a nontrivial solution of this system, i.e. corresponding determinant must be zero. 
This means that the following quadratic form is degenerative: 



M{r k Ak) = E ^(p) ■ MP) = E J ^ p ) ■ mp) 



(D.4) 



p.L- 



In this way we obtain the alternative equivalent system of J^{P) ( |2.9| ), used in sections 2.2 and 2.3, 
with three types co- vectors ipi (P) ; 02 (P) > 03 (P) ■ The transition from the original linear problem <f> J 
to the co-linear problem ip* <-» J* is important, because one can't define an evolution for 4> (but hence 
there is the ambiguity in the definition of the evolution for ip*). 
Introduce the translations for ip%.: 



r a ■ r k (p) = r k (aP) , n ■ r k (p) = r k {bP) 



(D.5) 



Then system ( ) can be written as 

J* = wi ■ ipl - Ui ■ "02 + K {T a + UiWi T&) ■ V3 = , 
J'2 = n ■ ipl + ku 2 w 2 T~ ■ + i w 2 — u 2 Tb) ■ 0| = , 

J3 = -u 3 • -01 + w 3 T^ 1 ■ lp2 + K (T a + U 3 W 3 ) = , 

thus the following 3x3 block form of the "monodromy matrix" D arises: 

/ Wi(n) -f/i(n) nT a + kU x {n)W x {n)T b \ 



(D.6) 



D*(n) = 



k k J7 2 (n) ^2 (n) T- 1 W 2 (n) - E/ 2 (rc) T 6 
V -U 3 {n) W 3 (n)T^ nT a + KU 3 {n)W 3 (n) ) 



(D.7) 



where block notations U k (n) = diagp u k (n, P), W k (n) = diagp w k (n, P). 

As an e xample o ne can consider the case N = 1. The determinant of D gives immediately all the 
invariants ( A. 6 A. 10 ) of single R. 
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